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ABSTRACT 



We present the first sub-arcsecond resolution two-dimensional stellar kinemat- 



ics and X-ray observations of the prototypical starburst galaxy NGC253 which 
define the position and nature of the galactic nucleus. These observations com- 
prise some of the best probes of the central 300 pc of NGC253, the nearest massive 
galaxy undergoing a powerful starburst, and will allow us to gain more insight 
into the nature of the centers of starburst galaxies. We get an estimate of the 
stellar kinematic center location corresponding to an area of r ~ 1.2" centered 
~ 0.7" southwest from the radio core, and historically presumed nucleus, TH2. 
Newly processed Chandra data reveal a central point-like hard X-ray source (X- 
1) lying ~ 0.4" southwest from the kinematic center. Very accurate alignment 
between radio, infrared and X-ray sources in the nuclear region shows that TH2, 
the IR photometric center and X-l are not associated with each other. As the 
kinematic center is consistent with the positions of TH2 and X-l, and both could 
be a manifestation of nuclear activity, we consider the two as possible galactic 
nucleus candidates. Although TH2 is the strongest compact radio source in the 
nuclear region, it does not have any infrared, optical or X-ray counterparts. If 



1 Based on observations at the European Southern Observatory VLT (075.A-0250). 
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the kinematic center is associated with this source, by analogy we suggest that 
the nucleus of NGC253 resembles our Galactic Center SgrA*. On the other hand, 
X-l is a heavily absorbed object (iV# = 7.5 x 10 23 cm -2 ) only detected at en- 
ergies > 2 keV (L 2 -iokeV ~ 10 40 ergs s _1 ). If X-l is instead associated with the 
kinematic center, the nucleus of NGC253 is compatible with an obscured low 
luminosity active galactic nucleus (AGN) or a spatially resolved super star clus- 
ter (SSC) brightening up in X-rays most probably due to young supernovae or 
supernova remnants, a situation also observed in the nuclear starburst of M82. 
If no SSC is associated with the kinematic center, we conclude that NGC253 is a 
galaxy in which a strong starburst and a weak AGN (either TH2 or X-l) coexist. 
Results from few other high resolution studies of nearby starburst galaxies (e.g. 
M82, NGC5253, NGC4945) indicate that the AGN in these systems, if present, is 
always in the low luminosity regime. This may indicate that the onset of nuclear 
activity in galaxies is closely related with the occurrence of star formation, and 
that we are witnessing the emergence or dissaperance of an AGN. 



Subject headings: galaxies: starburst - galaxies: individual (NGC 253) - galaxies: 
nuclei - galaxies: kinematics and dynamics- X-rays: galaxies - infrared: galaxies 



Introduction 



It is observationally well established that star formation in starburst galaxies manifests 
itself in the production of numerous young massive star clusters of extraordinary luminosity 
and compactness, normally referred to as super star clusters (SSCs). In fact, the four most 
luminous circumnuclear starburst systems (M82, M83, NGC2 53 and NG C4945) account for 
25% of the high-mass (> 8 M ) star formation within 10 Mpc ( JHeckmanl 1 1 99 81 ) . Nearby star- 
burst galaxies provide thus a unique laboratory for investigating the mechanisms which cause 
enhanced star formation to occur and be sustained, as well as physical conditions analogous 
to intensely star- forming galaxies identified at high-z (e.g. Shapley et al. 2003). An intrigu- 
ing issue in studies of starburst galaxies is the apparent absence of a well defined nuclear 
source which could be associated with the center of the galaxy. This situation i s present in 
spatially resolved sta rbursts in ne arby galaxies: e.g. M8 2 (IMatsumoto et aU200lf ); NGC4945 
(iMarconi et aJbood ): NGC5253 JSummers et al.ll2004h : NGC253 (this work). Massive black 
holes (BHs) and/or nuclear star clusters (NSCs) are thought to reside in the centers of most 



galaxies, with the first being mostly identified in massive ellipticals ( JRichstone et al.lll998 
Graham! 120081 ) and NSCs in irregular and late spirals (e.g. Boker et al. 2002; Seth et al. 



2010). Although tens of SSCs exist in the central region of starburst galaxies, the pres- 
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ence of a NSC in these systems is very poorly constrained. Since the central few hundred 
parsecs of starburst galaxies are heavily affected by large and nonuniform dust extinction, 
the determination of SSCs properties (age, colors, luminosities, etc.), and thus the detec- 
tion of a putative NSC, has been difficult and uncertain. The prese nce of massive BHs in 



galax y centers is normally inferred when nuclear activity is observed (jFilippenko fc Sargent 



19891 ). Active Galactic Nuclei (AGN) are usually characterized by prominent point sources 
in i^-band. Starburst galaxies do not show any IR point-like emission at the position of 
the assumed galactic nucleus, which in the majority of these systems has been associated to 
either a powerful compact radio source or a hard X-ray peak. However, in several cases it is 
not possible to reject/accept the AGN hypothesis on the basis of radio/X-ray emission peaks 
alone, particularly when the three photometric centers (radio, IR, X-ray) don't coincide. For 
the nearest, best studied starbursts, possible galactic center candidates include: young su- 
pernovae (SNe), supernova remnants (SNRs), intermediate- mass BHs, inactive supermassive 
BHs, low-luminosity AGN, etc. (see e.g. Wills et al. 1999 and Stevens et al. 1999 for the 
case of M82, Marconi et al. 2000 for NGC4945 and Fernandez-Ontiveros et al. 2009 for 
NGC253). 

As the center of a galaxy is likely to coincide with the bottom of the galaxy potential 
well, independent estimates of the kinematic center position from kinematic data can help 
to determine which photometric center, if any, is the most promising candidate for galactic 
nucleus. These important dynamical studies had not previously been widely performed 
because high S/N spectra over a full 2D field with high spatial resolution are required, and 
these are still scarce. In this paper we provide the first estimate of the kinematic center of 
the prototypical starburst galaxy NGC253 using sub-arcsecond resolution two-dimensional 
stellar kinematics. This is combined with data at radio, infrared and X-ray wavelengths to 
examine the nature of the galactic nucleus itself. These observations comprise some of the 
best probes of the inner region of NGC253 and will allow us to gain more insight into the 
nature of the centers of starburst galaxies. 

As the nearest massive galaxy experiencing a strong starburst (3.94 Mpc, 1" = 19 pc, 
Karachentsev et al. 2003), NGC253 is one of the most studied objects in the sky, yet many 
aspects of its nuclear region remain poorly understood. One of the main unresolved issues 
is the precise location and physical nature of the center of the galaxy. The galactic nu- 
cleus has historically been associated to a powerful compact non-thermal radio source (TH2 
in Ulvestad & An tonucci 1997) with characteristics similar to other synchrotron sources 
in active galaxies ( Turner fc Holll985l ): very small size (< 100 mas), high brightness tem- 



perature (T > 40000 K at 15 GHz), and a relatively flat spectral index between 8.4 — 23 
GHz(a = —0.3, S u oc u a ). Follow up studies accross the electromagnetic spectrum, each lim- 
ited by the spatial resolution available at that time, have matched this source with i) H2O 
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maser emission (jNakai et al.lll995l ). ii) different knots of I R emission (IForbes et al.l |2000| ; 
Galliano et al.ll2005l ). and iii) the hard X-ray emission peak (jWeaver et al.ll2002l ). suggesting 
that TH2 may be a low luminosity active galactic nucleus (LLAGN). However, it is now ap- 
parent from recent high-spatial resolution observations at different wavelengths, and using 
diverse observational techniques, that this picture may be incorrect. The evidences are the 
following: 



VLA and VLBA observations at 23 GHz (IBrunthaler et al.l 120091 ) have shown that 
the high velocity H 2 maser is not associated with TH2 but with another nearby 
continuum source (TH4 in Ulvest ad fc Antonucci 1997) , which has been classified as a 
SNR embedded in an HII region (ILenc fc Tingayll2006l ). This suggests that a maser- 
AGN connection is unlikely and that the maser emission is almost surely related to 
star formation. 



Very high spatial resolution VLBI observations at 2.3 and 23 GHz (ILenc fc Tingay 
20061 ; IBrunthaler et al.l 120091 ) did not detect TH2 on mas scales. As this source is 
unresolved and very bright in VLA 5-23 GHz continuum maps (beam sizes > 100 
mas), this questions the presence of a possible AGN. 

Very accurate alignment at parsec scales has been found between the near-IR and radio 
sources by means of Adaptive Optics V LT/NACO IR data with a compar able resolution 
to those in existing VLA radio maps (jFernandez-Ontiveros et al.ll2009l ). Surprisingly, 
this new alignment shows that TH2 has no infrared counterpart. The newly revealed 
IR nature of TH2 shows no evidence of a central point source that would indicate 
nuclear activity. However, the non-detection of IR emission does not rule out other 
sce narios such as the existence of an i nactive supermassive b lack hole (SMBH) like Sgr 
A* (jFernandez-Ontiveros et al.ll2009l ; IBrunthaler et al.ll2009l ). 



In addition, one needs to be cautious when associating the har d X-ray peak to TH2. 
It is well known that NGC253 harbors ~ 60 comp act radio sources (lUlvestad &: Antonucci 



19971 ) and ~ 37 near-infrared point-like sources ( Fernandez-Ontiveros et al.ll2009l ). This 
significant change in morphology with wavelength has led to different registrations in the 
literature (see e.g. Galliano et al. 2005 and references herein). Since Chandra data were 
processed with 0.5" pixel _1 resulting in ~ 1" resolution and have an astrometric accuracy 
of ~ 1" rms, the hard X-ray peak could be a ssociated with severa l other discrete sources in 
the field. This is demonstrated in Figure 2 of I Weaver et al.l (120021 ) where a variety of IR and 
radio sources are contained within t he size of the hard X -ray peak. Moreover, these authors 
followed the alignment proposed by IForbes et al.l (120001 ) which, in light of the current high 
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spatial resolution NACO images presented in lFernandez-Ontiveros et al.l (120091 ). seems to be 
incorrect. Thus, it is still a topic of active investigation whether TH2 harbors an AGN, or 
even if it corresponds to the galactic center, or not. 



Based on the alignment proposed by lFernandez-Ontiveros et al.l (120091 ) one could point 
to the infrared photometric center (or IR peak) as possible AGN candidate. The IR peak lies 
~ 3.3" SW from TH2. It is in fact the brightest source at infrared w avelengths (1—20 /im) and 
emits ~ 10% of the galaxy's bolometric flux (jMohan et al.l 120021 ). The spatial distribution 
of the 3.3 /im polycyclic arom atic hydrocarbons (PAH) feat ure-to-continuum ratio reaches a 
local minimum in this region (JTacconi-Garman et al.l 120051 ). This indicates that the source 
produces strong ionization, probably destroying the PAH molecules by photo-dissociation. 
Supporting this interpretation is the fact that multiband tracers of io nized gas such as 



Pa/3, Br7, [Hell and [Nell] are seen to peak strongly at thi s position (iBoker et al 



1998 



Engelbracht et al.l 11998 : Keto et al.l 119991 Forbes et al.l l2000l : iTacconi-Garmanl 120061 ). In 
addition, the IR peak has an individual counterpart at radio wavelengths (TH7 in Ulvestad 
& Antonucci 1997) which also presents a relatively flat spectrum between 8.4 — 23 GHz 
(a = —0.21, S v oc u a ), but it is ~ 4 times fainter than TH2. Moreover, this strong IR 
source is located at the vertex of the rather collimated ionised-gas region seen in NGC253 
(see Figure 8 of Tacconi-Garman et al. 2005). All these characteristics make this object 
unique among the other sources in the central region of the galaxy and could be explained 
either by invoking a very young SSC or an AGN. 

NGC253 has been the subject of many detailed kinematical studies of its gaseous 
components. Optical rotation curves of the ionized gas have shown that, i n addit i on to 



ordered rotation, large non-circular motions exi st within the central region ( lUlrichl 11978 
Munoz-Tunon et al. Hl993tlArnaboldi et al.lll9951 ). However, because of problems of obscura 



tion at optical wavelengths, results from Ha and [N II] observations are difficult to interpret. 
In additi on, HI and H a rotation curves are incomplete: either they cover the outer parts 
r > 100" (|PenceJ|l980j), or they are confined to the innermost regions where peculiar features 
and asymmetries in the ionized gas are dominant. Most of the optical and IR observations, 
aimed at studying the kinematics of the molecular gas in the nuclear region of NGC 253 
(e.g., Prada et al. 1996; Engelbracht et al. 1998), have relied on measuring the velocity gra- 
dients along chosen position angles passing through the nucleus (e.g., major axis and minor 
axis). Interferometric radio observations at millimeter and centimeter wavelengths measure 
two-dimensional velocity fields which provide a more complete picture of the kinematics 
i n the nuclear region. Several radio observations in molecular emission lines such as CO 



dMauersberger et al.lll996l; ICanzian. Mundv fc Scoyilla Il988t iDas et al.ll200ll ; |Paglione et al. 



200 4]), CS (Peng et al. 19961). H CN ( Paglione et al.lll9951 ). or in hydrogen recombination lines 



( lAnantharamaiah fc Gosslll9961 ) exist in the literature. These two-dimensional measurements 
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reveal a complex but systematic velocity pattern in the central region. In addition to solid- 
body rotation, the obs erved velocity fiel ds indicate motions which may be due to a barlike 



potential in the center (jPeng et al.lll996l ) or a kinematic subsystem which may be caused by 
a past merger event (Anantharamaiah & Goss 1996). Although most of these observations 
have high spectral resolution, they typically have poor spatial resolution > 2.5" and are not 
always useful for probing the dynamics very close to the nucleus. 

Because of this situation, different authors assumed different centers for their rotation 
curves/velocity maps: either the center is derived assuming complete symmetry of the outer 
parts in their fields of view, or it is coincident with the IR peak, or with the radio core 
TH2, or with the surface brightness peak. The highest resolution measurement of the two- 
dimensional velocity field available to date is that of Anantharamaiah & Goss (1996), who 
observed the H92a recombination line with a beam of 1.8" x 1.0" and a velocity resolution 
of 54 km s -1 . After trying several positions near the H92a emission peak, the optimum 
center of the velocity field was found to be at {a, 5) 20 oo = 00 h 47 m 33 s .06, -25°17'18".3 (la 
error ~ 0.3"), which is ~ 2.0" southwest of the H92a line peak which in turn is coincident 
with TH2. This means that the most accurate estimate of the kinematic center based on 
gas kinematics available to date lies closer to the IR peak (~ 1.3") than to TH2. 

However, kinematic studies of galaxies focused on the gaseous component (ionized 
and/or molecular gas) have the potential for revealing the dynamical forces that work in 
the gas closer to the galactic nucleus, but they cannot trace adequately the gravitational 
potential. The gas represents only a fraction of the total mass of the galaxy and could be 
driven by non-gravitational perturbations such as interactions, warps, bars, etc. Therefore, 
the kinematic parameters of any nuclear galactic disk (center, position angle and inclination) 
can be greatly influenced by these perturbations. Thus the velocity field of stars is needed 
to properly trace the gravitational potential. 

Integral-field spectroscopy is an ideal tool for studying objects with complex morpholo- 
gies and kinematics, such as the nuclei of starburst galaxies. This paper is the first of two 
describing the conclusions from an integral-field spectroscopic study with SINFONI on the 
VLT to clarify the structure, kinematics and emission mechanisms of the central 150 pc star- 
burst region of NGC253. In this part we use the CO stellar absorption features in K— band 
to determine the stellar kinematic center of the galaxy. The stellar CO features have the 
advantage that they trace a component of the galaxy nucleus that is not affected by non- 
gravitational processes such as shocks or winds. We also obtained archival Chandra data 
and re- processed them using 0.125" pixel -1 . The X-ray images are carefully compared with 
available high spatial resolution radio and near-IR data. This multiwavelength approach 
combined with the estimate of the kinematic center is essential to find the true nucleus 
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of NGC253. In a forthcoming publication (JMiiller-Sanchez et al.ll2010bl ) we investigate the 
physical properties and kinematics of the super star clusters and the molecular and ionized 
gas in the nuclear starburst. 



This paper is structured as follows. In Section 2, we describe the SINFONI and Chandra 
observations and the processing we performed on the data. The method used to extract the 
kinematics from the SINFONI spectral features is presented in this Section as well as the new 
X-ray images and spectra. In Section 3 we examine the stellar velocity field and obtain an 
estimate of the position of the kinematic center. We discuss the implications of our results 
for the nature of the galactic nucleus in Section 4. The overall conclusions of the study are 
outlined in Section 5. 



2. Data Analysis and Results 
2.1. Near-IR SINFONI data 



2.1.1. Observations and data reduction 



The SINFONI instrument on the VLT (IBonnet et all 12004 : Eisenhauer et al.ll2003J ) was 
used without the adaptive optics module, resulting in slice widths of 0.25", a pixel scale of 
0.125" pixel -1 , and a field of view of 8"x8". Due to the high-inclination of the galaxy (~ 78°), 
this field size is not enough to cover the diameter of the circumnuclear starburst regions under 
study, hence mosaicing was necessary. To cover the central 30" along the photometric major 
axis, 8 different pointings were performed and the raw data were obtained from the ESO 
archive. The individual frames were combined to form a mosaic that was used for analysis. 
The SINFONI mosaic has a dimension of 30" x 30" and is centered at the position of the 
presumed galactic nucleus, TH2. 

The SINFONI data presented here were taken in X-band (1.9-2.4 /im) on the night 
of 28 August 2005 in 18 individual frames of 300s resulting in 1.5h total integration time. 
In this night the atmospheric conditions were excellent (seeing < 0.7"), which allows us to 
obtain a resolution of 0.5" as measured from the FWHM of the smallest knot of emission in 
i^-band. The FWHM spectral resolution is 68 km s _1 (R/^4400). The da ta were reduced 
using the SINFONI custom reduction package SPRED (jAbuter et al.ll2005l ). This performs 
all the usual steps needed to reduce near infrared spectra, but with the additional routines 
for reconstructing the data cube. Flux calibration was performed using a standard star 
observed next in time to the science frames, which yielded in a 0.5" aperture centered at the 
near-IR emission peak a i^-band magnitude of 12.34. In addition, flux calibration was cross- 



checked with NACO data in a similar aperture ( [Fernandez- Ontiveros et al.ll2009 ). Agreement 
between the different data sources was consistent to 15%. 



2.1.2. Spectral fitting 

The shape of the CO absorption bandheads in a galaxy spectrum arises from the intrinsic 
spectral profile of the stars themselves convolved with a broadening profile which carries the 
information about the kinematics. The use of template stars allows one to separate these two 
components. We have extracted the 2D properties of the CO absorption features by fitting a 
Gaussian convolved with a template stellar spectrum to the continuum-subtracted spectral 
profile of the CO bandheads at each spatial pixel in the data cube. For the CO bandheads 
in NGC253, we chose a template stellar spectrum of a late spectral type star, HD49105 
(spectral type KO), observed with SINFONI in the same instrument configuration (pixel 
scale and grating). We performed a minimisation of the x 2 i n which the parameters of the 
Gaussian (amplitude, center and width in velocity space) were adjusted until the convolved 
profile best matched the data (x 2 /dof ~ 1). The uncertainties were boot-strapped using 
Monte Carlo techniques, assuming that the noise is uncorrelated and the intrinsic profile 
is well represented by a Gaussian ( jDavies et al.l 120071 ). This method allowed us to obtain 
uncertainties for the velocity map in the range of (5 — 15) km s^ 1 . 

Thanks to the high S/N of our data we could make use of the three overtone bands 
of 12 CO at 2.294 (2-0), 2.323 (3-1), and 2.353 (4-2) fjm. This is demonstrated in Fig. HJ 
where we have plotted the spectral region around the K— band CO bandheads of several 
knots of IR emission in the nuclear region. We have extracted stellar flux, velocity and 
di spersion maps. Howev e r, in th is paper we only analyze the stellar velocity field and refer 
to 



Miiller- Sanchez et al.l (J2010bl ) for a detailed study of the other maps. 



2.1.3. Distribution of the K—band continuum emission 

A flux map of .K"-band continuum emission in the central 570 x 570 pc of NGC253 is 
presented in Fig. [TJ The image is centered at the position of the radio core, and presumed 
galactic center, TH2 (jUlvestad fc Antonuccilll997l ). Regi stration of sources is based on th e 
alignment between radio and near-IR images proposed by lFernandez-Ontiveros et al.l ( 120091 ). 
Relevant IR sources in this study are labeled as "IR-n" following the enumeration used 
by these authors as well, except for the IR peak. Here we briefly summarize the main 
morphological features of the K— band continuum to put them in context for the upcoming 
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kinematic analysis and discussion. 

Our SINFON I K— band map agrees well with previous near-IR studies of slightly lower 
spatial resolution ( jSams et al.lll994l ; Engelbracht et al.lll998t Forbes et al.ll2000l) and with the 



high-r esolution Adaptive Optics VLT/NACO images presented in lFernandez-Ontiveros et al. 



(120091 ) . At these wavelengths (1.9-2.4 /^m) the continuum emission of the nuclear region 
consists of multiple compact sources distributed along the disc of the galaxy in an apparent 
ring-like structure superposed on a diffuse emission component. The spatial resolution of 
our data and the availability of mosaic images at H— and K— bands providing a larger field 
of view compared to standard SINFONI observations, allow us to identify several discrete 
sources in the central 300 x 300 pc of the galaxy. A qualitative and quantitative analysis of 
the intrinsic properties of the stell a r cluste rs found in the central 300 pc of NGC253 will be 
presented in iMuller-Sanchez et al.l (J2010bl ). 



As can be seen in Fig. HJ the IR photometric center (the point of maximum emission 
in the wavelength range 1-20 /^m) is offset by about ~ 3.3" (60 pc) from the radio core, 
TH2. While the IR peak has an individual counterpart at radio wavelengths (TH7 in Ul- 
vestad fc Antonucci 1997), there is no discrete IR source associated with the radio nucleus 
(JFernandez-Ontiveros et al.ll2009l ). As it was discussed in the Introduction, this situation 
has originated a profuse debate about the precise location of the center of the galaxy and 
the nature of these two photometric centers. Thanks to the availability of sub-arcsecond res- 
olution integral field data we are able to provide a third independent estimate of the center 
of the galaxy: the stellar kinematic center. This will be derived in Sect. 13.21 and compared 
with radio, IR and X-ray images of the nuclear region in Sect. 14.11 



2.2. X-ray Chandra Data 



2.2.1. Observations and Data Processing 



Level 2 event data from the ACIS instrument were extracted from Chandra archive 
(ObsID 3931) and observed in 2003 September 19. The data were reduced with the CIAO 4.Cc 
data analysis system and the Chandra Calibration Database (caldb 4.0.Qj). The exposure 



1 http://cda.harvard.edu/chaser/ 
2 See http://asc.harvard.edu/ciao 
3 see http://cxc.harvard.edu/caldb/ 
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time was processed to exclude background flares, using the task lc_CLEAN.sl|j in source-free 
sky regions of the same observation. The net exposure time after flare removal is 76 ksec. 
The nucleus has no significant pileup. 

Chandra data include information about the photon energies and positions that were 
used to obtain energy-filtered images, and to carry out sub-pixel resolution spatial analysis. 
Although the default pixel size of the Chandra/ ACIS detector is 0.492", smaller spatial scales 
are accessible as the image moves across the detector pixels during the telescope dither, 
therefore sampling pixel scales smaller than the default pixel of Chandra/ ACIS detector. 
It allows sub-pixel binning of the images. We made 0.125" pixel -1 images with FWHM 
resolution of ~ 0.5" in two bands: 0.5-2 keV (hereafter soft band) and 2-10 keV (hereafter 
hard band). The absolute astrometric accuracy of Chandra data is 1" rms. 

In addition to the high-spatial resolution analysis we applied smoothing techniques to 
detect the low contrast diffuse emission. We app lied the adapt i ve sm oothing CIAO tool 



CSMOOTH, based on the algorithm developed by lEbeling et all ( 120061 ) . csmooth is an 
adaptive smoothing tool for images containing multiscale complex structures, and it preserves 
the spatial signatures and the associated counts as well as significance estimates. A minimum 
and maximum significance S/N level of 3 and 4, and a scale maximum of 2 pixels were used. 

The central sources detected in the image analysis were extracted to make the spectral 
analysis using circular apertures centered in the sources and with small radii of r ~ 0.5". This 
radius was chosen as trade-off between number of counts and contamination from background 
emission. We then extracted the background from source-free apertures around the hard X- 
ray emission regions. We subtracted the background from each source and performed the 
spectral analysis in the background-free sources. Response and ancillary response files were 
created using the CIAO mkacisrmf and mkwarf tools. To be able to use the \ 2 as the fit 
statistics, the spectra were binned to give a minimum of 20 counts per bin before background 
subtraction using grppha task, included in FTOOLSci The analysis of the spectral counts 
was performed using the software package XSPEC (version 12.4. Qj, Arnaud et al. 1996). 



4 see http://cxc.harvard.edu/ciao/download/scripts/ 
5 http://heasarc. gsfc.nasa.gov/docs/software/ftools 
6 http://cxc. heasarc.gsfc.nasa.gov/docs/xanadu/xspec/ 
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2.2.2. X-ray imaging 

Fig. [2] shows the newly processed ~ 0.5" resolution ACIS images of the nuclear region 
(17" x 17") of NGC253 in two energy bands: soft (0.5-2 keV) and hard (2-10 keV) X-rays. 
The images are centered at the position of source IR-32, unambiguously identified at infrared, 
optical and X-ray spectral ranges (see Sect. 14.11) . 

The soft X-ray image shows several knots of emission embedded in a diffuse component 
which originate s in the nuclear starbur st and extends southwest in an apparent cone of ion- 



ization. Indeed, IStrickland et al.l (120001 ) showed that the soft X-ray emitting gas distribution 
is consistent with a conical outflow of hot plasma, that is well collimated in a direction nearly 
perpendicular to the plane of the galaxy and extremely si milar in morphology and ext ent to 



the previously observed limb-brightened Ha outflow cone (ILehnert fc Heckmanlll996l ). 



In the hard X-ray image the emission is more compact and less diffuse. Three individual 
point-like sources are clearly identified: two lying in the nuclear starburst (X-l and X-2) and 
one located ~ 6" SE from the nuclear region which corresponds to X-3 (Fig. [2]). We use the 
"X-n" nomenclature because it denotes easily the number of individual point-like hard X-ray 
sources found in the field (Fig. [2]). X-l and X-2 are separated by ~ 2.3" (44 pc). While X-2 
has a counterpart at soft X-rays, X-l is only bright at energies > 2 keV, suggesting that this 
source is highly absorbed. Extracting an even harder image in the 4-10 keV band shows 
that indeed X-l is the brightest point-like source at higher energies. This is consistent with 
the spectra of the three sources (see Sect. 12.2.31 and Fig. [3]). 



2.2.3. X-ray spectroscopy 

We performed spectral analysis for the three hard X-ray sources identified in the right 
panel of Fig. We investigated five baseline models for the three sources, two of them 
only with one component: abs (power-law) and abs(MEKAL); and the other three compos- 
ite models: abs(power-law)+abs(MEKAL), abs(power-law)+abs(power-law) and abs(power- 
law+MEKAL)+abs(power-law). The best fit was selected as the simplest model in which 
we can reproduce the observed spectrum based on x 2 statistics. We tested whether a 
more complex model is statis ticaly justified by using f-statistic (ftest within Xspec). See 



Gonzalez-Martin et al.l (120091 ) for a detailed description of the models and the fitting process. 
The spectra of the three sources are shown in Fig. [3] and the model parameters are listed in 
Table [U 
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3. Stellar kinematics 

3.1. Evidence for a rotating stellar disc 

We have extracted stellar kinematics in the nuclear region of the galaxy using the 
method described in Sect. I2.1.2[ allowing for the first time a two-dimensional study of the 
stellar motions on scales enclosing the star-forming nuclear ring. Fig. H] shows the velocity 
field of the stars in the central 30" x 30" of the galaxy. As can be seen in this Figure, the 
stellar velocity field is very regular, typical of rotational motions. The major kinematic axis 
exhibits a velocity gradient from -80 km s _1 in the NE to 80 km s _1 in the SW and is 
well aligned with the photometric major axis of the nuclear disc. In addition, it presents 
azimuthal symmetry and the zero-velocity axis (or kinematic minor axis) is perpendicular 
to the kinematic major axis, as expected in regions dominated by rotation. The observed 
regular pattern of the stellar velocity field permits us to use a simple axisymmetric disc model 
to characterize the stellar kinematics in the nuclear region, and particularly, to constrain the 
location of the kinematic center. The fitting method and the derived model parameters will 
be discussed in the following Section. 



3.2. Kinematic modeling: Defining the kinematic center of the galaxy 

If a velocity map exhibits ordered rotation, such as the stellar velocity field of NGC253, 
the velocity can be well represented by circular motion in a disc. The standard method 
for deprojecting galaxies from kinematic data is to assume that the galaxy is made up of a 
number of rotating coplanar annuli. The velocity field is fitted with the equation: 

V(r, 0) = V sys + V c (r) cos(0 - a) sini (1) 

where r is the radius of a circular ring in the plane of the galaxy (or the semimajor axis 
of the ellipse in the plane of the sky), V sys is the systemic velocity, V c is the ring circular 
velocity, i is the ring inclination, a is the position angle of t he projected major ax is, and <fi is 
the azimuthal angle. This is the cosine law approximation (jKrajnovic et al.ll2006 ). A crucial 



element of our investigation, which is intrinsic to this model, involves choice of the kinematic 
center. Most kinematic studies of galaxies do not consider the center as a variable, as they 
know a priori its location. For the purpose of this study, the stellar velocity field shall be 
fitted with a single disc model with four free geometric parameters: the center (xo,y ), the 
position angle a, and the inclination i. 

The parameters of the rotating disc can be well determ ined via the kinemetry method 



developed and described in detail by lKrajnovic et al.l ( 120061 ). Briefly, kinemetry is an exten- 
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sion of surface photometry to the higher-order moments of the velocity distribution. The 
procedure operates by first describing the data by a series of concentric ellipses of increasing 
major axis length, as defined by the system center, position angle, and inclination. Nor- 
mally, the location of the center is the only parameter that is determined prior to running 
a kinemetric analysis; the latter two parameters (position angle and inclination) can either 
be determined a priori and used as inputs or be measured functions of semi-major axis 
length. At each radius, a small number of harmonic terms in a Fourier expansion is needed 
to determine the best-fitting ellipse: 

N 

K(4>) = A + y^ A n sin(n0) + B n cos(n(p) (2) 

n=l 

where K is the best-fitting ellipse at a given radius, N is the maximum fit order used, and 
is the e ccentric anomaly, which i n the case of discs corresponds to the azimuthal angle. 



Following iTrachternach et al.l (120081 ). a decomposition of the velocity fields up to third order 
(i.e., iV = 3) is sufficient to capture most of the signal. If the motion is purely circular, 
the observed velocity along that ellipse will be described by a cosine law. In this ideal case, 
the Bi term should be the only coefficient present in the harmonic expansion. Power in the 
other coefficients represents either deviations from circular motion (bars, inflows/outflows, 
etc.), or artifacts induced by the noise or incorrect geometric parameters (kinematic center, 
position angle and inclination). 

Our fitting procedure therefore consists of iteratively changing the four geometric pa- 
rameters (xo, yo, a, i) until the values of A , Ai, A 2) A3, B 2 , and B 3 are zero or close to zero. 
In other words, we use the kinemetry software to construct the best fitting rotating disk 
model of the stellar velocity field and determine simultaneously the kinem atic center, posi 



tion a ngle and inclination. In contrast to oth er previous kinemetric studies (iKrajnovic et al. 



20071 ; iBoker et al.ll2008l ; IShapiro et al.ll2008l ). in which the center is known a priori and the 
position angle and inclination are determined in two steps, the nature of our problem sug- 
gests that we need to calculate the four parameters simultaneously, as the location of the 
kinematic center is thigthly connected with the position angle of the kinematic major axis. 
Note that this procedure affects the determination of non-circular motions and could lead to 
slightly different geometric parameters a and i as for the cases in which the centers are known 
a priori. The position of the kinematic center is to first order that position which minimizes 
the non-circular motions in a tilted-ring fit and adopting this center position would then 
result in non-circular motions slightly smaller than the ones derived for the cases of known 
kinematic centers. 

The regular pattern of the stellar velocities allows us to consider the entire velocity field 
at once and solve for a global kinematic center, position angle and inclination. In fact, we 
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do not consider the possibility of having a kinematic center which varies with radius as it 
does not have a physical meaning. In more detail we performed the following steps: 

1. We masked out pixels in the velocity field having < 10% of the peak flux in the K— band 
continuum and those presenting velocities greater than |150| km s _1 . 

2. We assumed initial values of (x , yo, a, i) for the stellar disc in units of pixels for (x , 
Do) and degrees for (a, i), and kept them constant through the kinemetric analysis. 
The choice of these units was dictated by practical reasons of easily identifying the 
kinematic center in our SINFONI maps. 

3. The observed velocity field V was sampled along a set of ellipses regularly spaced in 
the semimajor axis and with constant (xo, yo, ex, i) as defined above. 

4. The velocities of every ellipse were fitted with a complete kinemetric expansion. 

5. We determined the agreement between the observed velocity field V and the disc model 
V m (the 2D image reconstructed from only the B\ term of the kinemetric expansion) 
as: x 2 — J2j [(Vm,j ~ Vj)/AVj} 2 where N is the number of pixels included in the fit 
and AVj are the measured errors of the velocity at each pixel j. The residuals map 
{V m — V) reflects the combined powers in the higher coefficients, which may appear as 
asymmetries in the velocity field. As our stellar velocity field shows ordered rotation, 
the residuals are a measurement of the error in the estimation of the disc parameters. 



6. We then performed a multidimension al minimization of the function x 2 { x o, yo, PA 



i 



using the AMOEBA function in IDL flNelder & Meadlll965i l. 



7. We finally selected the best fitting values and the corresponding errors from the pa- 
rameter space of x 2 (xo, yo, ex, i). The values of (xo, yo, ex, i) that give the lowest 
X 2 , X 2 m in, were considered to be the best-fit values and the confidence intervals were 
determined by a range of x 2 values greater than x 2 m in, defined as A% 2 = x 2 ~ X 2 m in- 
If it is assumed that the observational errors are normally distributed, then A% 2 also 
follows a x 2 probability distributi on with the number of degrees of freedom equa l to 



the number of model parameters (JLampton et al.lll976t Ivan der Marel et al.lll998f ) 



The result of the above procedure provided the best fitting disc parameters acceptable 
to the 99.7% confidence level (corresponding to a 3a uncertainty): the center is located at 
pixel xq = 129lg , yo = 104^ in our data, a = 50 ± 11°, and i = 77l^ 3 degrees. For 
reference, TH2 is located at pixel x = 124 and y = 107 in our SINFONI maps, which implies 
that the kinematic center is offset by Ax = 0.6" and Ay = —0.4" from the radio core. Notice 
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that the position of TH2 is consistent with the kinematic center estimate within the errors. 
The estimated absolute position of the kinematic center is thus (a, 5)2000 = 00 h 47 m 33 s .14, 
— 25°17'17".52 with a mean 3a uncertainty of r ~ 1.2". The stellar velocity field used for the 
modeling is presented in the left panel of Fig. while the best fitting disc model is shown in 
the middle panel. The residuals between the data and the best fitting model are displayed 
in the right panel. It is apparent from the residuals map that the observed velocities can 
be remarkably well described by the adopted simple thin disc approximation. No obvious 
regular patterns are noticeable in the residuals map, suggesting that the features observed 
in the residuals are only due to noisy pixels in the data. Furthermore, this map shows values 
between [—15, 15] km s _1 consistent with the measured errors of the velocity map. The 
best fit a and i a re in agreement with those from t he photometric nuclear ring and those of 
the galactic disc ( IPencelll980l ; iPaglione et al.l 120041 ) . Regarding to the location of the stellar 
kinematic center, we found that it lies in between the two original candidates for the galactic 
center, TH2 and the IR peak. However, the estimated confidence interval of our kinematic 
modeling clearly rules out the IR peak as kinematic center (they are separated by ~ 2.6" 
distance from each other). A more detailed discussion on the candidates for galactic nucleus 
and their physical properties will be presented in Sect. HJ 



4. The Nucleus of NGC253 
4.1. Alignment between radio, IR and X-ray images of the nuclear region 



Fig. |H] summarizes the main results of our investigation. In this Figure, an overlay of 
the continuum emission at radio, IR and hard X-ray wavelengths in the central 10" x 10" 
(~ 200 x 200 pc) of the g alaxy is shown. The colour s cale c orresponds to the high resolution 
3.8/xm NA CO image fromlFernandez-Ontiveros et all ( 120091 ) . The 2-cm VLA A-configuration 
map from iTurner fc Hoi ( 19851 ) is shown in white contours and the new 2-10 keV Chandra 
image from this work is presented in blue contours. The relevant radio, infrared and X-ray 
sources are indicated in the Figure, and the stellar kinematic center within approximately the 
3<r errors is marked with a crossed circle of radius r = 1.2". We use the nomenclature from 



Ulvestad fc Antonuccil (119971 ) for the radio sources (except for source 5.73-39.5), and the 
IR and hard X-ray sources are labeled following the nomenclatures proposed in Sect. 12.1.31 
and Sect. 12.2.21 respectively. The alignment b etween radio and near-IR images was done as 
proposed by lFernandez-Ontiveros et all (120091 ). The key alignment, however, is between hard 
X-ray and IR images. The first registration of sources was obtained directly from Chandra 's 
absolute astrometry. The absolute astrometric accuracy of the nuclear pointing is ~ 1" rms. 
This astrometric approximation could be improved through the identification of common 
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sources. We noticed that two of the brightest IR point-like sources in the nuclear region, the 
IR peak and IR-32, have individual counterparts in the X-rays (see Fig. [6]). The presence 
of these reference objects allowed us to tie the X-ray and IR images and achieve ~ 0.5" 
relative position accuracy in searching for radio and IR counterparts to the X-ray sources. 
The alignment was further cross-checked using optical and X-ray images. A three-color 
composite H ST /Chandra image of the central 20" x 20" (400 x 400 pc) region of NGC253 
is shown in Fig. [7J Green corresponds to V^-band (WFPC2/F555W), red to soft X-ray and 
blue to hard X-ray emission. The image is centered at the position of source IR-32, which 
is an outstanding feature at IR, optical and X-ray wavelengths, and a reference object in 
the alignment. We found that the extended soft X-ray emission southeast of the nuclear 
starburst is well correlated with the l/-band emission in this region. Furthermore, the most 
prominent dark dust lane running along the northwestern part of the nuclear starburst is 
spatially coincident in optical and soft X-ray images (see Fig. [7J). 

As can be seen in Fig. |6l within the 3<r error circle of the stellar kinematic center, 
there are ~ 15 objects detected at differe nt wavelengths: 6 radio sources, 8 IR knots of 



emission interpreted as young star clusters (jFernandez-Ontiveros et al.ll2009l ). and the hard 



X-ray source X-l. Having ruled out the IR peak as the kinematic center, the two next likely 
candidates for galactic nucleus are the photometric centers at radio and X-ray wavelengths 
(TH2 and X-l, respectively), the two consistent with the kinematic center location. Surpris- 
ingly, TH2 and X-l are not each other's counterparts. X-l is located ~ 1.1" SW from TH2 
well beyond the error of ~ 0.5" in the registration of the images. Not detecting any point- 
like hard X-ray source at the position of TH2 questions strongly its AGN nature, but does 
not discard it as the galact i c nuc leus with characteristics similar to Sgr A*, as proposed by 



Fernandez-Ontiveros et all (120091 ). However, Chandra results clearly point to X-l as AGN 
candidate. As our analysis cannot conclude whether the stellar kinematic center corresponds 
to TH2 or X-l, we consider the two as possible galactic nucleus candidates and discuss the 
implications in the following subsections. 

As X-2 lies ~ 2.6" (50 pc) SW from the kinematic center and X-3 is located ~ 6" (120 
pc) away from the nuclear region (see Fig. [6]), the discussion of their properties is out of the 
scope of this paper. However, they are still presented in Fig. [3] and Tabled] for completeness. 



4.2. TH2 as candidate for galactic nucleus 

The newly processed Chandra observations of NGC253, combined with most of the cur- 
rent high-spatial resolution studies of this galaxy at different wavelengths, provide sufficient 
evidence against the AGN nature of TH2 (see Sect. [Q. In fact, this source is only detected at 
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radio wavelengths bet ween 2.3 and 23 GHz and using baselines corresponding to synthesized 
beam sizes > 30 mas (jSadler et al.lll995t iLenc &: Tingayl 120061 ). At frequenc ies lower tha n 
2.3 GHz the source is not visible probably due to strong free-free absorption (JTingayl 12004 ) . 
Longer baselines completely resolve TH2 as demonstrated by recent VLBI observations at 
2.3 and 23 GHz with 15 x 11 and 2.5 x 1.1 mas resolution, respectively (ILenc fc Tingay 
20061 ; iBrunthaler et al.ll2009l ). The upper limit on the flux density of these two observations 
is ~ 1.1 mJy. As this source is unresolved with the VLA (size < 100 mas and a flux density 
of 20 mJy at 23 GHz) , this implies that its linear size is probably in the range of 0.5 — 2 pc. 
Lenc &: Tingayl ( 120061 ) performed LBA observations at 2.3 GHz with slightly shorter base- 
lines detecting TH2 with a flux density of ~ 6 mJy, a source size of 80 x 60 mas (1.6 x 1.2 
pc), and a brightness temperature of ~ 10 5 K. The spectrum of the source between 5-15 
GHz is relatively flat (a ~ —0.47) and can be best fit with a free-free absorbed-power-law. 
All these characteristics are not particular unusual in the radio regime and can be explained 
in terms of either an acc reting massive object similar to other synchrotron sources in the 



centers of active galax ies (IHummel et al. 



198. 



a very compact SNR (jUlvestad fc Antonucci 



19971 ). 



Turner fc Holll985l : iGallimore et al.ll2004l ). or 



The apparent lack of detected source variability at rad io wavelengths over 8 year tim escales 
implies that if TH2 is a SNR i t is likely not expanding (JUlvestad fc Antonuccilll997| ). Re- 
cently, IBrunthaler et al.l (120091 ) reported signs of possible variability of TH2. However, as 
pointed out by these authors, a comparison of different values of flux densities in the liter- 
ature is difficult, first due to potential contamination by the extended diffuse emission from 
the galaxy, and second because of the different resolutions used in the observations. This 
uncertainty ca n be seen in the comparison o f the flux densities from VLA A and B configu- 
ration data in lUlvestad fc Antonucci! (119971 ) and in Table 1 and Figure 5 of ILenc fc Tingay 
( 120061 ) . The different flux density measurements for TH2 have originated confusio n on the 
nature of its radio spectrum, which has been argued to be flat (ITurner fc Holll985l ). 



Brunthaler et al.l ( 120091 ) question the presence of a possible LLAGN based on the non- 
detection of TH2 on mas scales with an upper limit of 1 mJy at 23 GHz. In their analysis, they 
suggest that TH2 may be a Gigahertz Peaked Spectrum (GPS) source due to an apparent 
turnover in its spectrum between 4.8 and 8 GHz. However, as they pointed out, the GPS 
nature of TH2 is dubious as GPS sources having a turnover point between 5-10 GHz normally 
have sizes around 20 pc, ~ 10 times larger than the upper limit of 2 pc predicted for TH2, 
and are several orders of magnitudes more powerful than TH2. In addition, the lack o f any 
optical or IR counterpart is inconsistent with the GPS hypothesis ( jSnellen et al.lll998l ). 



In summary, TH2 is the brightest radio source (a factor > 3 compared to other sources) 
detected by the VLA (beam sizes > 100 mas) in the central 150 x 150 pc of NGC253. 
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However, it does not have any infrared, optical or X-ray counterpart, which argues a gainst 
an AGN or SNR nature. Regarding the last possibility, Fernandez- Ontiveros et all (120091 ) 
show that a scaled version of the Crab or N40 SNRs to TH2 radio emission level should 
produce an IR/optical detectable counterpart, which is not the case. If indeed associated 
with the kinematic center of NGC253, these properties make TH2 similar to the center of 
our Galaxy, SgrA*. A dormant BH like SgrA* at the distance of NGC253 (3.95 Mpc) would 
have flux densities at radio (23 GHz) and IR (flares at 3.8 //m) wavelengths of ~ 3.6 /xJy 
and ~ 0.02 /iJy, respectively, making the sour ce undetectable. The existence o f a scaled 
version of SgrA* in NGC253 was proposed by IFernandez-Ontiveros et al.l (120091 ) based on 
the IR/radio ratio of TH2, measured with VLA 15 GHz and NACO L— band images of 
comparable resolution (~ 200 mas, 4 pc), which resembles that of SgrA* rather than an 
AGN. 



4.3. X-l as candidate for galactic nucleus 



Chandra observations reveal a heavily absorbed source of hard X-rays embedded within 
the 2>a error c i rcle o f the kinematic center (see Fig. [6]). The hard X-ray emission map from 
Weaver et al.l ( 120021 ) shows a partly extended component along the nuclear ridge confined to 
an encircled area of r~ 2.5". At these scales, it is not clear whether or not the peak in hard 
X-rays is pointlike, but it has been argued that it is coincident with TH2. We reprocessed 
Chandra data and found, with an accurate astrometry of ~ 0.5" rms, a central hard X-ray 
peak (X-l) that is indeed point-like but which lies ~ 1.1" SW from TH2. It is evident from 
the similarity of the morphologies between the two hard X-ray maps that X-l corresponds 
to the strongest source in the 4.5-8 keV continuum map of IWeaver et al.l ( 120021 ). 



While the X-r ay images from the t wo studies are very consistent, the X-ray spectrum 
of X-l obtained by I Weaver et al.l ( 120021 ) is very different to the one measured in this work 
(see Figure 3 of Weaver et al. 2002 and our Fig. [3]). We believe that this discrepancy is 
caused by the different aperture sizes used in the extraction of the spectra. In their work, 
a nuclear spectrum was extracted from within a circular region of diameter 5" centered 1" 
to the northeast of TH2 in order to avoid contamination from X-2. However, this aperture 
is with high probability contaminated by the diffuse 0.5-8 keV X-ray emission surrounding 
X-l (Fig. ED- 

The X-ray spectrum of X-l presented in Fig. [3] is best fitted with two absorbed power- 
laws with N H1 = 4 x 10 22 and N H2 = 7.5 x 10 23 cm" 2 (see Sect. EJL3] and Table UJ). While 
Nhi is consisten t with the extinction of the nuc lear starburst region measured with near- 
IR observations (IFernandez-Ontiveros et al.l 120091 ). iY# 2 is particularly high, indicating that 
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the source is heavily absorbed. This is supported by the lack of emission in the soft X- 
ray band. We hardly detected photons with energies below 2 keV resulting in a very low 
F .5-2keV / F2-iokeV ratio of 3 x 10~ 3 . X-ray s ources with low F()K_ 9 . kpV / F 9 ._ mkf y and high 



N H are most often regarded as hidden AGN (jUeda et al. 1 120071 ; IWinter et al.ll2009l ). We 



estimate the intrinsic 2- 10 keV luminosity o f X-l to be 1 x 10 ergs s , consistent with the 



lower limit proposed by (IWeaver et al.ll2002l ) of 2 x 10 39 ergs s _1 . The intensity of the X-ray 
emission places the source in the "ultraluminous X-ray source" category. Given its location, 
coincident with the kinematic center and embedded in the nuclear starburst, chances are 
high that the ultraluminous source is a central massive object, either a weakly accreting 
massive BH (a LLAGN, Terashima et al. 2002; Gonzalez-Martin et al 2009) or an evolved 
NSC brightening up in X-rays most probably d ue to young SNe and /or SNRs, similar to the 
luminous X-ray source near the center of M82 (IStevens et al.lll999l ). We examine these two 
possibilities below. 



Within the achieved 0.5" alignment error between radio, IR and X-ray images (Sect. H~Tj) . 
there are two SSCs that could be associated with X-l: IR-11 and IR-18 (Fig. [6]). The two 
SSCs are spatially resolved in L— band, and are very similar in siz e (FWHM ~ 3 — 5 pc), flux 
densi ty (~ 2 — 2.5 mJy in L-band) and shape of their IR SEDs (jFernandez-Ontiveros et al. 



20091 ). While IR-11 has an individual counterpart at radio wavelengths (TH6 in Ulvestad 
& Antonucci 1997), there is no discrete radio source associated with IR-18 (see Fig. [6]). 
In the radio regime, TH6 is the second strongest flat-spectrum source in NGC253 (a ~ 
—0.01, Lenc & Tingay 2006). Notice that it is not clear whether the spectrum of TH2 is 
flat or not (see Sect. 14. 2p . It appears to be resolved in VLA observations at 23 and 15 
GHz (sizes of 4.5 x 2.4 and 4 x 1.8 pc, respectively), but unresolved at 8.4, 5 and 2.3 
GHz. It was not de tected at 1.4 GHz with the LBA probably due to free-free absorption as 
TH2 (jTingayl 12004] ). Interestingly, this source shows a position offset in LBA observations 
at 2.3 G Hz that is not as con sistent as the other sources when compared against VLA 
positions ( JLenc fc Tingayl 120061 ) . These facts suggest that TH6 may have structure that 
varies significantly with frequency. In addition, it is associated with Radio Recombination 
Lines (RRL) emission that is consistent with the free-free absorbed spectra of HII regions. As 
i n young; stellar clusters, hard X-ray radiation is main ly produced by young SNe and/or SNRs 
( IStevens et al.lll999t I Van Bever fc Vanbeverenl 120001 ). the SSC IR-ll/radio-TH6 is the most 
likely counterpart of X-l, presenting the expected simultan eous contribution in radio, IR an d 
X-ray from young SNe and/or SNRs evolving in the SSC (JFernandez-Ontiveros et al.ll2009l ). 
The fact that IR-18 has no radio counterpart makes more diffi cult to justify it s assoc iation 
with X-l. There is an additional radio source, 5.73-39.5 in iLenc fc Tingayl ( 120061 ). also 
compatible with the location of X-l (Fig. E]). However, its flux density has been measured 
only at two frequencies: 8.4 and 23 GHz. In these measurements 5.73-39.5 is four times 
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weaker than TH6, and has no IR counterpart, all together making dubious any interpretation. 

A second alternative is the association of X-l with a hidden LLAGN on the basis of the 
hardness of its X-ray spectrum with an absorbing column density of almost 10 24 cm -2 . No 
other manifestation of the active nucleus is found, neither a strong radio nor a near-infrared 
point source is located at the position of X-l. The absence of any signature other than in the 
hard X-rays implies that the AGN is non-sta ndard and similar to the hidden AGN detected 



in NGC4945, but ~ 100 times less luminous (Marconi et all 120001 ). 



In summary, X-l appears to be a heavily absorbed object (Nh > 10 23 cm -2 ) only 
detected at energies > 2 keV (L 2 _iofceV ~ 10 40 ergs s _1 ). If indeed associated with the 
kinematic center of the galaxy, the simplest explanation accounting for these properties is 
that X-l is a hidden LLAGN, similar to the one observed in NGC4945. Alternatively, within 
the Chandra error circle of r ~ 0.5", X-l could be associated with two SSCs (IR-18 and IR- 
ll/radio-TH6) and one resolved radio source (5.73-39.5). Among these three possibilities, 
the SSC IR-ll/radio-TH6 is the most likely counterpart to X-l. In this case, the hard X-ray 
emission would be associated to processes occuring in the SSC which could produce the 
intrinsic 2-10 keV luminosity of 10 40 ergs s _1 , being the most plausible explanation young 
SNe and/or SNRs evolving in this massive cluster. The nature of the kinematic center of 
NGC253 would the n be a NSC of 10 5 — 10 6 Mm with properties similar to other clusters in 



the nuclear region ([Fernandez- Ontiveros et al.ll2009l ; iMiiller- Sanchez et al.ll2010bf ). 



X-l might alternatively be a background X-ray source. In this case, the inferred lu- 
minosity is only a lower limit and a background AGN could well be invoked to explain the 
Xray power of the source. Although we cannot completely rule out this scenario, it seems 
very unlikely that there exists a background AGN which is precisely coincident with the 
kinematic center of NGC253. In addition, as it was discussed above, X-l may have radio 
and/or IR counterparts which are believed to belong to the galaxy. 



5. Conclusions 

SINFONI near-infrared integral field spectroscopy combined with Chandra imaging and 
spectroscopy have been used in this work to probe the galactic nucleus of NGC253. Thanks 
to our high S/N integral field data we were able to extract 2D stellar kinematics in the 
nuclear region of NGC253 and estimate the kinematic center of the galaxy with a very good 
accuracy. Our main results and conclusions can be summarized as follows: 

• The stellar kinematic center is located in between the strongest compact radio source 
(TH2) and the IR photometric center, or IR peak. It lies ~ 0.7" SW from TH2 and 
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~ 2.6" NE from the IR peak. 

• The estimated confidence interval of our kinematic modeling (r ~ 1.2") clearly discards 
the IR peak as kinematic center, but includes the position of TH2. 

• Newly processed Chandra images reveal a central point-like hard X-ray source (X-l) 
lying ~ 0.4" SW from the kinematic center. 

• Very accurate alignment between radio, IR and X-ray sources in the nuclear region 
shows that TH2 and X-l are not associated with each other, being separated by ~ 1.1". 
As the two could be a manifestation of nuclear activity, and our current estimate of 
the kinematic center location is consistent with both sources, we consider the two as 
possible galactic nucleus candidates. 

• TH2 is the brightest radio source, with the highest brightness temperature, detected 
with the VLA in the central region of NGC253. However, it does not have any IR, 
optical or X-ray counterparts. When observed at very high spatial resolution with the 
VLBI, TH2 is among the weakest sources (~ 6 mJy at 2.3 GHz, Lenc & Tingay 2006), 
or it is not detected (< 1 mJy at 23 GHZ, Brunthaler et al. 2009). If the kinematic 
center is associated with TH2, we conclude that the nature of the nucleus of NGC253 
resembles that of SgrA*. 

• X-l is a heavily absorbed object (Nh > 10 23 cm~ 2 ) only detected at energies > 2 keV 
(L 2 _iofceV ~ 10 40 ergs s _1 ). If X-l corresponds to the kinematic center of the galaxy, 
the simplest explanation accounting for these properties is that the nucleus of NGC253 
harbors a hidden LLAGN. 

• Alternatively, X-l may be associated with one SSC in this region: SSC IR-11/radio- 
TH6. In this case the kinematic center of NGC253 would then be a 10 5 — 10 6 M 
spatially resolved (diameter of ~ 5 pc) NSC in which some of its massive stars are 
evolving into the SN phase. 

An important conclusion from this study is that the nature of dusty galaxies cannot be 
always characterized by photometric data alone but must be complemented with spectro- 
scopic observations. In particular, we demonstrated the power of integral field spectroscopy 
for the determination of kinematic centers of galaxies. Further near-IR integral-field obser- 
vations with high spatial resoultion would be useful to derive a more accurate position of 
the kinematic center, measure directly the mass of the putative central massive object, and 
reveal the presence of either broad recombination lines or forbidden high ionization lines. A 
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detection of any of these would provide definite evidence of the existence of an AGN in this 
galaxy. 
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Fig. 1. — SINFONI K— band continuum emission map of the central 30" x 30" (570 x 570 
pc) of NGC253. The image is centered at the position of the radio core, and presumed 
galactic nucleus, TH2 (black cross). This corresponds to an absolute position of (a, 5) 2 ooo 
= 00 h 47 m 33M7, -25°17'17".l. Figures 1, 4 and 5 (all SINFONI maps) are centered at this 



positio n. Relevant IR sources are labeled as "IR-n" following Fernandez- Ontiveros et al. 



( 120091 ) enumeration, except for the IR peak. The spectra of the five IR sources indicated in 
the Figure in the wavelength range between 2.27 and 2.4 /im are shown. The three overtone 
bands of 12 CO at 2.294 (2-0), 2.323 (3-1), and 2.353 (4-2) /im, used for the extraction of the 
stellar kinematics are clearly seen. The spectra were extracted within a circular aperture of 
r = 0.25". 
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Fig. 2. — Chandra soft (0.5-2 keV, left) and hard (2-10 keV, right) X-ray images of the 
central 17" x 17" region of NGC253. In the two panels north is up and east is left. The 
images are centered at the position of source IR-32, unambiguously identified at infrared, 
optical and X-ray wavelengths. This corresponds to an absolute position of (a, 5)2000 = 
00 h 47 m 33 s .4, — 25°17'14".0. Figures 2, 6 and 7 (all Chandra images) are centered at this 
position. The locations of the three hard X-ray sources discussed in the paper are marked 
in the right panel following the nomenclature "X-n" . 
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Fig. 3. — Spectral fits to the hard X-ray sources X-l (top), X-2 (middle), and X-3 (bottom) 
See Table 1 for model parameters. 
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Fig. 4. — Two dimensional velocity field of the stars in the central 30" x 30" of NGC253. 
The velocity map is centered at the position of TH2 (black cross) and the black contours 
delinate the K— band continuum emission. North is up and East is left. 
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Fig. 5. — Maps showing the data and the best-fit kinematic model used to determine the 
kinematic center of the galaxy. The maps are centered at the position of TH2 (black cross) 
and the contours delinate the K— band continuum emission. Left: Stellar velocity field used 
for the modeling, Middle: Best-fit kinematic model of a rotating disc, Right: Residuals 
(data-model) to the fit. 
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Fig. 6. — Nuclear region of NGC253 at infrared NACO L-band (3.8 /im, colour scale), 
radio VLA-A 15 GHz (white contours), and hard X-ray 2-10 keV (blue contours) wave- 
lengths. Relevant radi o, infrared and X-ray sources are indicated in the Figure. We use 
the nomenclature from lUlvestad fc Antonuccil (119971 ) for the radio sources "THn" (except 
for source 5.73-39.5). Infrared s ources are labeled as "IR-n" following the enumeration of 
Fernandez-Ontiveros et al.l (120091 ) (except for the IR peak), and the X-ray sources are marked 
as "X-n" . The 3<r confidence interval of the stellar kinematic center location is approximated 
by a crossed circle of radius r = 1.2". 
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Fig. 7. — Three-color composite image of Chandra soft X-ray (red), hard X-ray (blue) and 
HST V— band (green) emission in the central 20" x 20" region of NGC253. The image is 
centered at the position of source IR-32, which is an outstanding feature at IR, optical and 
X-ray wavelengths, and a reference object in the alignment. 
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Table 1: Spectral models for the hard X-ray sources found in the nucleus of NGC253. 
Source Type a Counts b N m c N H2 C T F , 5 _ 2keV d F 2 _ wkeV d L 2 _ 10keV c X 2 / dof 



X-l 

X-2 
X-3 


2po 

po 

po 


393 

436 
803 


4 0+ 3 ' 7 

9 cc+0.8 

7 +0 - 3 

u -' -0.2 


7 cr+24 
'°-12 


3.6±?:| 

9 q+0.4 
z -°-0.3 

1 fi +0 - 2 


-15.0+?i 

-i4.2±8:l 

-13.718:1 


-12 64+ - 03 

i o i 7+O.O8 
J- - 1 ' -0.35 

-12 Q8+ - 05 


40.03+?:° 2 

38.18_ ; 33 
38.2918:? 6 


15.0/15 
33.1/18 
32.3/43 























a Bcst-fit spectral model, po stands for abs(power-law), me for abs(MEKAL), and 2po for abs(power- 
law)+abs(power-law). 

6 Number of counts in the range 0.5-10 keV 
c Hydrogcn column densities in units of 10 22 cm~ 2 
^Observed flux in logarithmic units of ergs cm~ 2 s _1 
e Intrinsic hard X-ray luminosity in logarithmic units of ergs s _1 (assuming d = 3.95 Mpc) 

Note. — All quoted untertainties are 90% except for the fluxes and luminosities, which are la. 



